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Abstract 


We  show  by  temperature-programmed  desorption  that  hydrogen  desorbs  from  Ge(100)2xl 
near  570  K  with  near-first-order  kinetics,  similar  to  the  behavior  of  hydrogen  on  Si(100)2xl.  Hie 
near-first-order  desorption  kinetics  are  attributed  to  pairing  on  surface  dimers  induced  by  the  it  bond 
on  unoccupied  dimer  atoms,  and  a  pairing  enthalpy  of  5±1  keal/mol  is  inferred.  However,  a 
comparison  between  the  pairing  enthalpies  for  H  atoms  on  Ge(100)  and  Si(100)  with  the  electronic 
structure  of  the  respective  clean  surfaces  indicates  that  estimates  of  the  it  bond  strength  based  on  the 


surface  band  structure  do  not  correlate  with  the  propensity  for  pairing. 
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The  kinetics  and  dynamics  of  hydrogen  desorption  from  Si(100)2xl  have  been  a  major  focus 
of  recent  effort  Hydrogen  chemisorption  has  proved  a  useful  complement  to  studies  of  geometric 
and  electronic  structure  of  semiconductor  surfaces  and  plays  a  critical  role  in  the  growth  of  Si  and 
GexSii-x  thin  films  by  chemical  vapor  deposition  and  atomic  layer  epitaxy.  Much  of  the  recent  work 
was  stimulated  by  the  discovery  that  recombinative  hydrogen  desorption  follows  first-order  kinetics 
on  Si(lOO),1*2  in  contrast  to  the  second-order  behavior  seen  on  metal  surfaces.  The  original  propo¬ 
sal  by  Sinniah  et  al.1  that  the  first-order  kinetics  results  from  the  rate-limiting  excitation  of  hydrogen 
atoms  to  a  band-like  delocalized  state  was  contradicted  by  observations  of  near-second-order  desorp¬ 
tion  kinetics  of  hydrogen  from  Si(l  1 l),2-4  by  measurements  of  the  diffusion  kinetics  of  hydrogen 
on  Si(lll)  which  suggest  conventional  hopping,5  and  by  dynamical  measurements  showing  that 
desorbing  H2  is  rotationally  cold  and  vibrationally  hot,  implying  a  highly  symmetric  transition 
state.6  Wise  et  al.2  suggested  that  the  first-order  desorption  kinetics  are  due  instead  to  pairing  of 
hydrogen  on  the  dimerized  surface  atoms,  and  Boland7  and  we8  proposed  that  preferential  pairing  of 
H  atoms  is  a  consequence  of  the  it  bond9*10  on  “unoccupied”  dimers.  Boland  obtained  direct 
evidence  for  preferential  pairing  by  scanning  tunneling  microscopy  (STM)  and  estimated  a  it  bond 
strength  of =18  kcal/mol  from  tunneling  spectra  above  un-  and  singly-occupied  dimers.7  We 
reviewed  existing  evidence  for  it  bonding  and  preferential  pairing  and  proposed  a  doubly-occupied 
dimer  model  to  quantitatively  describe  the  desorption  kinetics.8  The  driving  force  for  pairing, 

A//pair,  is  equal  to  the  difference  in  Si-H  bond  strengths  between  doubly-occupied  dimers,  with  two 
hydrogen  atoms,  and  singly-occupied  dimers,  with  one  hydrogen  atom,  and  was  identified  with  die 
it  bond  strength8  by  analogy  to  molecular  it  bonds.11  Atfpair  was  estimated  as  -7.5  kcal/mol  by 
comparing  model  predictions  to  temperature-programmed  desorption  (TPD)  data1*2  indicating  slight 
departures  from  first-order  kinetics  9  HOfer  et  al.  demonstrated  that  the  kinetics  do  indeed  deviate 
from  first-order  at  low  coverage,  and  used  our  model  to  obtain  a  refined  determination  of  6±1 
kcal/mol  for  AHp^.12  A  modest  value  for  Affpair  is  supported  by  high-level  ab  initio  calculations  of 
Si-H  bond  strengths  in  a  cluster  analogue  of  Si(100)2xl  that  imply  values  of  2-5  kcal/mol.13*14 

The  interaction  of  hydrogen  with  Ge(100)2xl,  which  has  received  much  less  attention, 
provides  a  test  for  the  predicted  generality  of  preferential  pairing7*8  and  also  offers  insight  into  Ac 
relationship  between  AF/pair.  the  it  bond  strength,  and  the  surface  electronic  structure.  The  behavior 
of  hydrogen  on  germanium  is  qualitatively  very  similar  to  that  on  silicon.  At  coverages  up  to  one 
monolayer  a  (2xl)  H  monohydride  structure  is  formed;15*17  a  less-stable  dihydride  can  be  formed  at 
higher  coverages.17  H2  desorbs  from  Ge(l  1 1)  with  near-second-order  kinetics,18  just  as  on 
Si(l  1 1).2*4  Adsorption  of  hydrogen  on  Ge(100)  was  investigated  in  early  work18  but  TPD  results 
were  not  reported  We  found  that  H2  resulting  from  decomposition  of  H2S  or  H2O  on  Ge(100) 
desorbs  with  approximately  first-order  kinetics  near  570  K,  implying  an  activation  energy  for 
desorption  near  40  kcal/mol.19  The  rate-limiting  step  in  [100]  growth  of  germanium  from  Gc2H^, 
presumably  H2  desorption  from  Ge(  100)2x1,  was  found  by  reflectometry  to  follow  first-order 
kinetics  with  an  activation  energy  of  40  kcal/mol.20  In  the  present  work  we  demonstrate  that  the 


desorption  kinetics  of  H2  from  Ge(100)2xl  show  systematic  deviations  from  first-order  behavior  in 
the  low  coverage  limit  which  can  be  quantitatively  described  by  the  doubly-occupied  diner  model. 
Comparison  of  the  inferred  value  of  A//pair  with  the  surface  electronic  structure  of  Ge(100)  shows 
that  estimates  of  the  n  bond  strength  based  on  the  band  structure  of  clean  Ge(100)  and  Si(100)  do 
not  correlate  with  the  propensity  for  pairing  of  H  atoms.  However,  experiments  on  nearly- 
hydrogen-saturated  surfaces  might  allow  for  a  spectroscopic  determination  of  the  x  bond  strength. 

Experiments  were  performed  in  an  ultrahigh  vacuum  (UHV)  chamber  equipped  with 
LEFD/ESDIAD  optics,  an  Auger  spectrometer,  a  quadrupole  mass  spectrometer  (QMS)  with  a 
water-cooled  shroud,  an  ion  gun,  and  a  tungsten  filament  for  dosing  of  atomic  hydrogen.19  The 
Ge(100)  sample  was  cleaved  from  a  0.25-mm-thick  n-type  wafer,  p  =  5-40  £2  cm,  to  a  dimension 
of  13.4  mm  x  13.8  mm.  A  chromel-alumel  thermocouple  was  cemented  into  a  small  hole  drilled 
near  one  edge.  After  degreasing,  the  Ge  sample  was  placed  in  the  chamber  and  cleaned  by  several 
sputter-and-anneal  cycles  (i'at*  =  2-4  |xA  cm-2,  £at+  -  500  V,  Tm*,* 1  -  850  K). 

Atomic  hydrogen  exposures  were  performed  by  backfilling  the  UHV  chamber  with  H2  to 
pressures  of  2xlO'8-5xlO'7  Torr  and  heating  a  coiled  W  filament  located  “3  cm  from  the  sample  to 
1700-1800  K.  After  dosing,  when  the  background  pressure  fell  to  2-3  x  10-10  Torr,  the  sample 
was  rotated  to  face  the  entrance  slit  of  the  QMS  shroud,  at  a  distance  of  »  0.5  cm.  The  temperature 
was  ramped  at  2  K  s*1  and  the  QMS  signal  for  m/e=2  was  recorded  by  an  AT-compatible  personal 
computer.  A  coverage  calibration  for  hydrogen  was  obtained  by  TPD  of  H2S,  whose  adsorption  as 
H  +  SH21  saturates  at  0.5  monolayer,22  and  which  yields  exclusively  H2  and  GeS  upon  heating.19 

TPD  traces  for  H2  following  atomic  hydrogen  exposures  are  shown  in  Fig.  1  for  various  initial 
coverages,  6b.  The  peak  desorption  temperature,  Tp,  at  all  initial  coverages  is  about  570  K.  The 
nearly  coverage-independent  values  of  Tp,  together  with  the  asymmetric  peak  shape  at  high  6b, 
indicate  near-first-order  desorption  kinetics.23  However,  the  slight  but  readily  discernable  increase 
in  Tp  and  the  more  symmetric  peak  shape  at  Iowa  6b  demonstrate  a  departure  from  first-order 
kinetics.8  The  dependence  of  Tp  on  6b  is  shown  in  Fig.  2  together  with  model  predictions. 

The  doubly-occupied  dimer  model  comprises  an  exact  solution  to  a  lattice  gas  model  which 
incorporates  preferential  pairing  of  adsorbates  on  dimers  but  neglects  interactions  between  dimers.8 
The  coverage  of  hydrogen  which  is  paired  on  surface  dimers  is  given  by8*24 


©2=6- 


_ 26(1-6) _ 

(l  +  4[exp(AHpai^K7>l]©(l-6)}1/2  +  1 


(1) 


where  6  is  the  instantaneous  total  coverage  of  hydrogen,  R  is  the  gas  constant,  and  T  is  the 
temperature.  The  assumption  of  quasiequilibrium,  i.e.,  that  hydrogen  diffusion  is  fast  compared  to 
desorption,  is  supported  by  measurements  of  diffusion  and  desorption  kinetics  of  hydrogen  on 
Si(l  1 l)4*5  and  Si(100).12  At  sufficiently  low  coverage  diffusion  may  become  rate-limiting,  but  tins 


feature  is  omitted  from  the  mode'.  Recombinative  desorption  is  assumed  to  occur  between  H  atoms 
paired  on  surface  dimers, 2,6-8,12  so  that 

2  (2) 

Assuming  that  the  preexponential  factor  for  k  is  2xl015  sec*1,  as  was  found  for  H2  desoiption  from 
Si(100)  by  Hofer  et  al .,12  who  took  H-atom  pairing  into  account  in  the  analysis,  the  TPD  peak 
temperature  at  high  initial  coverage  implies  an  activation  energy  of  42  kcal/mol.  The  dependence  of 
Tp  on  60  predicted  by  the  model  is  shown  in  Fig.  2  for  several  assumed  values  of  AZ/pair.  By 
comparing  the  experimental  and  model  results,  we  infer  AZ/pair  to  be  5±1  kcal/mol  for  H  on 
Ge(100),  which  is  only  slightly  less  than  the  6±1  kcal/mol  obtained8*12  for  H  on  Si(100).  The 
inferred  value  of  AHpah  is  quite  insensitive  to  the  assumed  value  of  the  preexponential  factor.  If  the 
preexponential  factor  is  assumed  to  be  1013  sec*1,  the  activation  energy  must  be  reduced  to  36 
kcal/mol  to  be  consistent  with  the  value  of  Tp  at  Gq  =1.  The  values  of  Tp  then  predicted  for  AZ/pa^ 
=  5  kcal/mol  are  only  slightly  higher  than  those  shown  in  Fig.  2  and  are  in  equally  good  agreement 
with  the  data.  At  6b  =  0. 1,  for  example,  the  predicted  value  of  Tp  increases  by  only  0.4  K  by 
changing  the  preexponential  factor,  which  is  well  within  the  error  bars  of  »2  K  associated  with  the 
uncertainty  of  ±1  kcal/mol  in  AHp^. 

STM  has  not  yet  been  applied  to  hydrogen  adsorption  on  Ge(100),  but  infrared  spectroscopic 
measurements16  provide  independent  evidence  for  at  least  a  degree  of  preferential  pairing  of 
hydrogen.  Infrared  absorption  features  at  1979  and  1991  cm*1  for  Ge(100)2xl:H  at  one-monolayer 
hydrogen  coverage  arise  from  the  asymmetric  and  symmetric  stretch  modes  of  hydrogen  atoms  on 
doubly-occupied  dimers  and  display  distinct  polarization  properties.16  Distinct  s-  and  p-polarized 
spectra,  with  a  two-peak  structure  in  the  latter,  were  also  observed  for  nominal  exposures  as  low  as 
10  L,16  for  which  we  estimate  a  hydrogen  coverage  of  0.3  ML.  If  occupation  of  the  dangling  bond 
sites  were  random,  then  rally  30%  of  the  hydrogen  would  have  been  present  as  doubly-occupied 
dimers  and  the  spectrum  would  have  been  dominated  by  a  one-peak  structure.8  A  degree  of 
preferential  pairing  could  occur  during  adsorption  at  a  temperature  where  the  adsorbed  atoms  are 
immobile  by  either  a  “hot  precursor”  mechanism25  or  by  abstraction.  If  localization  of  adsorbing  H 
atoms  requires  several  bounces  to  dissipate  the  energy  released  by  formation  of  the  bond  to  the 
surface,  adsorbing  atoms  will  tend  to  settle  preferentially  into  the  deepest  wells,  viz.,  on  already- 
singly-occupied  dimers.  Abstraction  of  surface  hydrogen  and  halogen  atoms  by  incident  atomic 
hydrogen  appears  to  be  quite  facile,  at  least  on  Si(lOO).1*26  If  abstraction  occurs  even  while 
adsorption  is  taking  place  at  low  coverage,  abstraction  of  more  weakly-adsorbed  hydrogen  (at 
singly-occupied  dimers)  will  occur  preferentially  to  that  of  more  strongly-adsorbed  hydrogen  (at 
doubly-occupied  dimers)  due  to  a  smaller  activation  barrier,  which  would  also  lead  to  an 
enhancement  in  the  population  of  doubly-occupied  dimers. 


The  similarity  in  the  values  of  A/Zpaj,-  f°r  hydrogen  on  Ge(100)  and  on  Si(100)  is  consistent 
with  the  similarity  in  it  bond  strengths  in  H2Ge=GeH2  and  H2Si=SiH2  (=25  kcal/mol).27  The 
chemical  similarity  of  germanium  and  silicon  also  extends  to  the  geometric  and  electronic  structure  of 
their  (100)  surfaces.  The  basic  structural  units  on  reconstructed  Ge(100)  are  buckled  dimers,28'31 
with  several  ordered  arrangements  possible.  The  dimer  bond  length  may  be  estimated  as  2.41  A 
from  the  measured32  value  of  its  parallel  component  (2.34  A)  and  the  calculated31  buckling  angle 
(14*).33  This  value  is  intermediate  between  the  bond  length  in  molecular  digermenes,  2.21-2.35 
A,34  and  the  single  bond  length  of  2.44  A  in  bulk  germanium,  and  indicates  a  highly  strained  double 
bond.  Similar  considerations  hold  for  Si(100),  although  the  question  of  buckled  versus  symmetric 
dimers  has  been  controversial.  The  preponderance  of  recent  work35'39  indicates  that  the  dimers  are 
instantaneously  predominantly  buckled.  The  Si=Si  dimer  bond  length,  measured35  as  2.32±0.02  A 
and  calculated13*40*41  as  2.19-2.40  A,  is  similarly  indicative  of  a  highly  strained  double  bond.  STM 
images  of  buckled  dimers  on  Ge(lOO)28  and  Si(lOO)42  appear  quite  similar:  the  amplitudes  of 
occupied  and  unoccupied  states  are  greatest  on  the  “up”  and  “down”  atoms,  respectively,  of  the 
dimer  pair.  The  dimer  tilt  angle  is  significantly  greater  on  Ge(100)  (=14*)  than  on  Si(100), 
measured35  as  7*  and  calculated40  as  7-10*.  Both  the  larger  tilt  angle  for  germanium  dimers  and  the 
apparently  larger  barrier  to  flipping28*31*40  are  analogous  to  the  deviation  from  and  barrier  to 
planarity  in  trans-bent  H2Ge=<JeH2  versus  H2Si=SiH2«43 

What  is  the  7C  bond  strength  in  surface  dimers  on  Ge(100)  and  Si(100)?  Perhaps  the  cleanest 
definition  for  it  bond  strengths  in  molecules  is  die  activation  energy  for  rotation  about  die  double 
bond,1 1*27  which  is  inapplicable  to  surface  dimers.  A  second  measure  is  based  on  the  heat  of 
hydrogenation  of  a  double  bond1 1  and  is  equivalent  to  our  definition  of  AZ/pan-.8  For  H2C=CH2 
and  H2Si=SiH2  the  two  definitions  of  the  it  bond  strength  agree  to  within  a  few  kcal/mol.11  A  less 
direct  definition  is  the  energy  difference  between  the  ground  state  singlet  and  the  lowest-energy 
triplet  configurations,  which  also  agrees  to  within  a  few  kcal/mol  with  the  rotation  barriers  for 
H2Si=SiH2  and  H2Ge=GeH2.11<27  However,  the  band  structure  of  dangling  bond  surface  states 
reflects  the  periodicity  of  die  surface  and  provides  no  obvious  measure  of  the  it  bond  strength  or  of 
the  propensity  of  H  atoms  to  pain  bond  formation  is  intrinsically  a  local  process  and  chemisorption 
disrupts  the  surface  periodicity  and  band  structure.  A  simple  localized-it-bond  picture  of  the  clean- 
surface  electronic  structure  is  contradicted  by  the  dispersion  of  both  the  occupied  (it  or  Dup)  and 
unoccupied  (it*  or  Dfown)  bands  in  the  T-T  direction  (parallel  to  the  dimer  rows),  with  bandwidth  s 
of  0.8  and  0.7  eV  on  Si( 100)2x1 3639,44  and  1.0  and  0.35  eV  on  Ge( 100)2x1 ,293°  respectively. 
The  average  band  gap  over  the  surface  Brillouin  zone  is  significantly  larger  for  Ge(100),  =1.9 
eV,30*45  than  for  Si(100),  *  1  eV,39.46  showing  the  opposite  trend  to  that  of  A/fpair  and  of 
molecular  it  bond  strengths.  Boland  estimated  the  it  bond  strength  in  a  slightly  different  way,  as 
twice  the  difference  in  occupied-dangling-bond  binding  energies  of  unoccupied  and  singly-occupied 
dimers,7  the  analogue  of  which  gives  a  reasonable  estimate  of  the  it  bond  strength  in  H2Si=SiH2.47 
However,  the  dangling-bond  electronic  states  probed  by  the  STM  are  likely  to  be  strongly 
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delocalized,  as  STM  spectral  measurements  of  the  surface  band  gap  yield  similar  results  to  those 
measured  by  conventional  means30*39*46  on  both  Si(lOO)7*10  and  Ge(lOO),28  and  are  therefore  of 
questionable  applicability  for  determining  a  “local”  7t  bond  strength.  While  a  delocalized  interaction 
could  in  principle  supply  the  driving  force  for  pairing,  a  simple  count  of  the  singly-  and  doubly- 
occupied  dimers  in  Boland’s  STM  image  of  H  on  Si(100)  after  annealing  to  630  K7  leads  to  an 
estimate  of  6-7  kcal/mol  for  A//pair  using  Eq.  (I),48  entirely  consistent  with  the  value  inferred  from 
the  desorption  kinetics8*12  but  inconsistent  with  Boland’s  higher  estimate7  for  the  7t  bond  strength. 

We  conclude  that  the  estimate  of  the  it  bond  strength  represented  by  AH pair,  5±1  kcal/mol  for 
H  atoms  on  Ge(  100)2x1  and  6±1  kcal/mol  for  hydrogen  on  Si(lOO),9*13  quantitatively  accounts  for 
preferential  pairing  of  hydrogen  and  that  a  reliable  it  bond  strength  cannot  be  readily  extracted  from 
the  band  structure  of  the  clean  surface.  However,  it  would  be  intriguing  to  investigate  the  electronic 
structure  of  isolated  un-  and  singly-occupied  dimers,  surrounded  by  doubly  occupied  dimers  and 
prepared  either  by  careful  dosing  of  atomic  hydrogen  and/or  by  brief  annealing  of  a  hydrogen- 
saturated  surface,7  to  see  whether  a  spectroscopic  measure  of  the  it  bond  strength  could  be  obtained. 
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FIG.  1.  Temperature-programmed  desorption  (TPD)  data  for  H2  following  exposures  to  atomic 
hydrogen.  The  initial  hydrogen  coverage  in  monolayers,  60,  is  indicated  for  each  trace. 
The  vertical  line  indicates  the  TPD  peak  temperature  for  0q  =  1. 
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FIG.  2.  Dependence  of  TPD  peak  temperature  of  H2  on  initial  hydrogen  coverage.  Experimental 
results  are  shown  together  with  predictions  of  the  doubly-occupied  dimer  model  as  a 
parametric  function  of  AH^T.  Tp  was  calculated  using  an  activation  energy  of  42  keal/mol 
and  a  preexponential  factor  of  2xl015  s’1  for  k  and  a  heating  rate  of  2  K  s'1. 


